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The launch vehicle dynamics affected by bending and sloshing modes are considered. 
Attitude measurement data that are corrupted by flexible modes could yield instability of the 
vehicle dynamics. Flexible body and sloshing modes are reconstructed by sliding mode 
observers. The resultant estimates are used to remove the undesirable dynamics from the 
measurements, and the direct effects of sloshing and bending modes on the launch vehicle 
are compensated by means of a controller that is designed without taking the bending and 
sloshing modes into account. A linearized mathematical model of Ares I launch vehicle was 
derived based on FRACTAL, a linear model developed by NASA/MSFC. The compensated 
vehicle dynamics with a simple PID controller were studied for the launch vehicle model that 
included two bending modes, two slosh modes and actuator dynamics. A simulation study 
demonstrated stable and accurate performance of the flight control system with the 
augmented simple PID controller without the use of traditional linear bending filters. 
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Nomenclature 

pitch angle 

pitch angle rate 

translational velocity 

the i th bending mode (z = 1 , . . . , n) 

the rate of the i th bending mode 

the j th sloshing mode (j = 1 

the rate of the j th sloshing mode 

the state variables of the actuator 

the state vector of the actuator 

the control function (the input of the actuator) 

the output of the actuator (the input of the vehicle model) 
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I. Introduction 

D ue to significant structural mass constraints such that the payload mass ratio is maximized, launch vehicles tend 
to be highly flexible. In-flight winds cause bending moments that negatively affect vehicle strength, while 
bending motion coupled with sloshing of propellants produces stability problems. Much research in 1960s was 
directed at the control of launch vehicles for stable, high-performance atmospheric flight 1,2 . Traditional control 
designs have accommodated highly flexible launch vehicles, including the Saturn V 2 3 ; while more recently, new 

approaches using theory have been applied for the Ariane launch vehicle. 4 Attitude control of a flexible launch 

vehicle using adaptive notch filtering in a cascade with real-time PID control has also been addressed. 5 The 
mentioned techniques are based on filtering the flexible mode components via shaping a frequency characteristic of 
the launch vehicle 1 ' 4 or notch filtering the excited oscillations. 5 However, filtering for attenuation or phase 
stabilization of flex modes can add phase lag near the control frequency, affecting rigid body performance, and 
design of classical bending filters for robustness to large bending mode uncertainties is a significant challenge. 

Completely new approaches based on estimating bending and sloshing modes with consequent “cleaning up” 
the corrupted measurements and compensating the coupling terms were proposed for the first time in the work by 
Shtessel and Baev. 13 The proposed technique used a novel sliding mode observer that operates in parallel with a 
traditional PID controller. The first bending mode was taken into account while designing the flight control system. 
The proposed 13 technique is extended in the present paper, the contribution of which can be outlined as follows: (1) 
two different sliding mode observation techniques are studied and used for observing bending and sloshing modes, 
(2) two bending modes and two sloshing modes are taken into account, and (3) the proposed techniques are applied 
to a coupled linear model of Ares I crew launch vehicle (CLV) that includes attitude dynamics, two bending and 
sloshing modes, and an actuator model. Accurate performance of the proposed Ares I attitude control system that is 
robust to both two bending and sloshing modes as well as to the actuator dynamics has been confirmed via computer 
simulations. 


II. Concept 

The idea of the proposed method is in estimating the bending and sloshing modes and using the estimate to 
generate the control. The proposed algorithm comprises the following steps: 

(1) The measurements corrupted by bending mode deflections and rates are to be cleaned up so that primarily 
rigid body measurements remain, and the clean measurements are to be used in the controller. 

(2) The coupling terms from bending and sloshing modes that affect the CLV dynamics are to be compensated 
via the control functions that use the coupling term estimates. 

As a result, the flexible CLV compensated dynamics are expected to behave as rigid body compensated dynamics. 

Two different sliding mode observers are proposed. They are 

(1) A sliding mode disturbance observer. In this case the bending modes are treated as disturbances, with 
unknown dynamics. The advantage of this approach is in its robustness to the bending mode dynamics. 

(2) A sliding mode state observer. In this case all states, including bending and sloshing ones, are observed 
together, and then the necessary compensations are performed. The advantage of this approach is the capability to 
simultaneously clean up the measurements and compensate for the coupling. The disadvantage is in requiring the 
exact knowledge of the CLV dynamic model including the bending and sloshing modes. 


III. Fundamentals of sliding mode observers 


A. Classical state sliding mode observer 16 

Consider linear time invariant completely observable system 
[ x = Ax + Bu 
\y = Cx 


( 1 ) 


where x GU n , y (=U m ,u (=U m , A(=U nxn , B GU nxm , n>m . The output y is assumed measurable. The goal is to 
estimate x so that as time increases. System (1) is transformed by a nonsingular transformation 


Tx = 


xi 

y 


( 2 ) 
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to a regular form 


1*1 =A n x l +A n y + B l u 

\y = ^ 21 * 1 + ^n.y + ^2 U 


(3) 


The following sliding mode observer (SMO) is proposed; 

I Xy — A j + A Y2 y “1“ + Tv 

y = A2iXi+A 22 y + 52 M - v 


(4) 


where v s § is the injection term that is to be designed, and LeJ <n m)xm is the matrix that is to be determined in 
order to provide e x - x x - x x — » 0, e y = y - y — » 0 . 

The estimation error dynamics are derived as 


[e 1 - A 1*1 + \l e y + ^ v 

[*y = A 21 e l +A 22 e y ~ v 


(5) 


and the sliding mode injection term 

Vf = Msign(yi - ), M > 0, i -1,2,..., m (6) 

drives system (5) to the sliding mode ( ^ = 0 ) 


f*l “(At + LA 2l) e l 

K=° 

in finite time. The matrix L is to be determined in order to provide desirable eigenvalues for the compensated error 
dynamics (7). Definitely, such matrix exists, since the system (1) is observable. 


Remark 1. The sliding mode observer (4) is equivalent to the Luenberger observer of reduced order. 


B. State and Disturbance Sliding Mode Observer 6 

Consider a linear, time invariant, completely observable system 
(x = Ax + Bu + Dfi(x,t) 

\y = Cx 

where x eD ”, y e □ ”, u e □ m , A e □ nxn , B e □ nxm , n > m , and /j (xj) e L /c is an unknown bounded disturbance. 

The goal is to estimate x so that as time increases, and reconstruct the disturbance f\(x,t) . 

System (8) is transformed by a nonsingular transformation (2) to a regular form 
| A =A n x l +A l2 y + B l u 
\y = A 21 x 1 +A 22y + B2 u + D 2 f\ 


It is assumed that the matrix A n 41 ^ is Hurwitz. 

Consider the SMO of the form 

1 x x — A 1 ^ x i + A^ 2 y B\U — A^2&y 

y - At A + A 22y + b 2 u + v 

The estimation error dynamics are derived as 

| A = A i*i 

[ e y = A 2\ e \ + A 22 e y ~B>2f\ +v 


where e y = y-y, e\ = X] -x\ are the estimation errors. 
The following sliding mode injection term 


(z? + L)6 


p > 0 


drives e y — » 0 in finite time with |a 2 i^ + A 2 2 e y -D 2 /i < L . 


( 10 ) 


( 11 ) 


( 12 ) 
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Indeed, by introducing a Lyapunov’s function candidate for second equation of eq. (1 1) 


and requesting 


V = Q.5e{,e v 


VZ-p\\e y 


we obtain by substituting the injection term (12) into the derivative of the Lyapunov’s function candidate: 

T t / \ || || y || || T \P + 

v - e yZy - e y \^2\ e \ + A 22 e y ~ ^2 A + V j - | e y + V = | e y 1 | — = 


Remark 2. Fulfilling inequality (14) yields a finite convergence time e y = y - y — » 0 . 

Therefore, the estimation error dynamics becomes in finite time 
|A - \ i e i 

n ( 16 ) 

b=° 

and ei —> 0 as time increases. 

The disturbance term D 2 /l can easily reconstructed as time increases using a concept of equivalent control 
v eq =D 2 fi (17) 

where v eq is to be estimated by low-pass filtering of the high frequency switching injection term 

v eq =LPF(v) (18) 

C. State and Measurement Noise Sliding Mode Observer 6 

Consider the linear time invariant completely observable system with the measurement corrupted by the bounded 

noise f 2 (x,t) gD m 

f x = Ax + Bu 

^ , (19) 

[y = Cx + f 2 (x,t) 

where v g □ n , y g □ m , u g □ m , A g □ nxn , B # fj^ xm , n > m . The goal is to estimate x so that x^x as time 
increases, and reconstruct the noise of measurement f 2 (x,t) . 

System (19) is transformed by a nonsingular transformation (2) to a regular form 

1*1 = A 1*1 + A ny~ A2/2 + B i u 
{ y - Ai*i + A22 y + b 2 u + f 2 

It is also assumed that the matrix A n g □ i s Hurwitz. 

Consider the SMO of the form (10). The estimation error dynamics are derived 
\h = A n e \ + A nf 2 

\ . • ( 21 ) 

\fy “ ^21^1 + A 22 e y + A + v 

The sliding mode injection term (12) drives e y —> 0 in finite time with |a 21 ^ +A 22 e y +/ 2 || < L . The observation 
error dynamics become, in the sliding mode 

= A n e \ + A2/2 

\ • ( 22 ) 

[A - ~ A 2l e l ~ v eq 

Excluding Cj from eq. (22) we obtain a dynamic filter for f 2 reconstruction 


f 2 (s) = - s! + A 2l (SI-A ll ) 1 A l2 v eq (s) 
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D. State observer and disturbance reconstruction: general case 


12 


Consider the nonlinear, minimum phase 
(pit) ©0 m 

= f(x) + G(x)(p(t) 

I }’ = h(x) 


completely observable system with the bounded disturbance 


(24) 


where xeD n ,yeD m ,/(x)<=D n , G(x ) ©}J nxm ,n>m . The goal is to estimate x so that as time increases, 

and reconstruct the bounded disturbance <p(t) . The state observation and the disturbance reconstruction technique, 
based on differential geometry 15 and higher order sliding mode control (HOSM), 7 is developed for system (24) in 
Ref. 14. 


IV. Linearized Mathematical Model of Flexible Crew Launch Vehicle 


The linear model of the CLV used herein is derived from a subset of the linear stability analysis package 
FRACTAL 17 (Frequency Response Analysis and Comparison Tool Assuming Linearity), developed for control 
design and sensitivity studies at NASA/MSFC. The FRACTAL model includes the effects of bending mode 
coupling with propellant slosh, interaction due to the motion of attached elements such as rocket nozzles, and a 
comprehensive aerodynamics model. 

The model originates from a Lagrangian approach to the system dynamics which results in a set of coupled 
linear equations of motion describing vehicle motion with respect to a pseudo -inertial reference frame, with the 
assumption of small perturbations about an operating condition. The fundamental assumptions are that the vehicle 
has tetragonal inertia symmetry, that it ascends on a gravity turn trajectory, that the roll rate is constrained to zero, 
and that mass properties are quasi-constant over the analysis interval, and that structural flexibility can be well- 
represented via the superposition of a set of orthogonalized flexible modes with the rigid-body dynamics. 

The derived model used for the present study is presented in a form of LTI system which accounts for pitch 
angle, pitch rate, translational velocity, n bending modes with corresponding rates, and m sloshing modes with 
corresponding states. The model provides two outputs, a pitch angle (e) and pitch rate ( q ), the measurements of 
which are corrupted by the bending mode dynamics. 

The model can be described by the LTI system 
fx = Ax + BU 


y -Cx 


(25) 


and the actuator dynamics are also presented in a form of an LTI system 
f z = Q\Z + 0. 2 x + Pu 
[U = Q. 3 z 


where x = (d,q,v,rj l ,rj 2 ,...,T ?2n _ l ,r ?2n ,S l ,S 2 ,...,S 2m _ l ,S 2m ) T eD 2 " +2m+3 , 

Since both coordinates of interest (i.e. 0,q) are available for measurements, the output vector y gQ can be 
presented as a superposition of these two coordinates and all aforementioned bending modes, such that 

y — Cjjq + C 2 x 2 (27) 

where Xi =(0,qf «fjl 2 , x 2 =(v,rj l , S h S 2 ,...,S 2m _ l ,S 2m f e □ 2n+2m+1 . 


Remark 3 . It is worth noting that Cj is the identity matrix, and the first and last 2m + 5 columns of the matrix C 2 
are (0,0)^ . 

The control input u is a scalar output of a PID-like controller that is to be designed so that 

lim | >’ic 0) - >i 0)| = lim 1 0 C (r) - d(t)\ = 0 . (28) 

t — >00 t — >00 


V. Flexible CLV attitude control using sliding mode disturbance observer 

In this section we will design the sliding mode disturbance observer based control upon the following 
assumptions: 
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Al. Only one bending and one sloshing mode are taken into account 
A2. The dynamics of actuator are neglected. 

The mathematical model (25) can be partitioned and rewritten as follows: 

(a) Vehicle attitude dynamics 

0 - q 

< q — Cf 2 \0 + $23 V + ($24/71 a 25 r l2 ^ 26^1 ® 21^2 ) 4" ^21^1 

V — $21^ 4- $33V + ($34/7^ + $35/72 + $3^<5| + $37^2 ) 4" ^31 w l 

(b) Bending mode dynamics 

|7l=% 

1% =fl[ 51^ + a 53 v + a 54 7 7l + c, 55 ? 72 + a 56^'l +a 57^2 +^51 M 1 + +»! 

(c) Sloshing mode dynamics 

^1=^2 

^2 ~ $71^ “1“ $73V + $74/71 + $75/72 4" $76^1 4" $77^2 ^71^1 

(d) Measurements 

=6» + q 4 7i 

3^ =<? + $25 7 72 

The problem can be formulated as 

(A) Given flexible CLV dynamical equations, design a simple PID controller that provides 6{t) ^ 6 c {t) as 
time increases without taking into account coupling terms from slosh and bending modes as well as 
neglecting bending mode effects on attitude measurements. 

(B) Design a Higher Order Sliding Mode (HOSM) disturbance observer that estimates bending modes and 
“cleans up” the measurements. 

(C) Improve the control system robustness to the perturbations by retaining the perturbed system stability and 
restoring the stability margins. 


(29) 

(30) 

( 31 ) 

(32) 


A. PID controller design 


The PID controller is designed for the simplified attitude dynamics given by eq. (29) assuming 
that v = tji = q 2 ~ = $2 = ^ • The ITAE-based controller that drives eg = 0 C — 0 — >0 as time increases in 

accordance with differential equation 

2*3 

Cq +1.75 co^cq +2.15 co n CQ + co n eq =0 (33) 

is presented 


u \ =tH4 

b 2 \ v 


+ \.15co n 0 c 


- a 2 \0 -l. 1 5 co n q + 2. 15 co^eQ + $^{ eqdt^j 


(34) 


It is worth noting that in order to implement the PID control law in eq. (34) the commanded pitch rate and 
acceleration terms 0 C , 0 C must be computed. It can be done either numerically or using HOSM differentiators 7 . 

B. Simulations of the compensated flight control system of Ares I without cleaning up the measurements 

Simulation results of the system (29), (32), (34) with perfect measurements y@ =6, y q = q and no coupling 
between eq. (34) and eqs. (30), (31) with co n = 1 are shown in Figs. 1 and 2. 

The results of the simulations of the coupled system (29)-(32), (34) with perfect measurements yq =0, y q = q 
and a> n = 1 are shown in Figs. 3 and 4. 
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Fig. 1 Ideal output tracking 
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Fig 2 Ideal control function 
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Fig. 3a Output tracking in system with coupling and 
perfect measurements 



Fig. 3 a Zoomed output tracking in system with coupling 
and perfect measurements 



Fig 4a Control function in system with coupling and 
perfect measurements 



Fig 4b Zoomed control function in system with coupling 
and perfect measurements 
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C. Simulations of the compensated flight control system assuming exact rate measurements, i.e. y q = q . 

The results of the simulations of the coupled system (29) -(3 2), (34) that takes into account bending and sloshing 
mode effects with perfect rate measurement y q = q and corrupted attitude measurement while co n = 1 are shown in 

Figs. 5 and 6. The simulations of the coupled system (29)-(32), (34) that takes into account bending and sloshing 
mode effects with corrupted both attitude and rate measurements demonstrated instability. 



.010 

.008 




Fig. 5a Output tracking in system with coupling, exact rate 
measurement and corrupted attitude measurement 


Fig 5b Zoomed output tracking in system with 
coupling with coupling, exact rate measurement and 
corrupted attitude measurement 



Time (sec) 

Fig 6a The control function for system with coupling, exact 
rate measurement and corrupted attitude measurement 


Fig 6b Zoomed control function for system with 
coupling exact rate measurement and corrupted 
attitude measurement 


D. Bending mode observer design using sliding mode observation technique 

In this section the bending mode in eq. (30) is treated as unknown disturbance for the CLV dynamics (29). 
Equation (29) is rewritten in a new basis 
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( 35 ) 


y e =0 + c u rj x 


8 =ye-cu r h 

< y q = ^ + ^25 7 72 

-> < 

q = y q -C25 V2 and < 

V = V 

V 


V - V 


8 = ye- c um 
q = yq-c 2 5 1 ii 

V = V 


This is 


ye = y q +{ p 1471-^572) 

‘ y q — a 2 \ y$ + (~ a 2\ c \4 +a 24 ) 7 1\ + a 23 v + a 25 T l2 + ^25 , h + a 26 $l + a 21^2 + W\ ll \ (36) 

v = a 2l y G +a 33 V + (-a 31 c 14 + 034)771 + a 25 r] 2 +a 26 8 x + a 2l^2 + ^ 3 i M i 
The state observer for system (36) is designed as 

ye = y q +v \ 

< y q =« 21)'0 + a 23 v + O l (.) + v 2 (37) 

V = +^33V + 0 2 (.) 

where the terms ®,(.)=a 26 ^ + a 21 S 2 +b 21 u l , 0 2 (.)=a 36 ^ + a 37 <S 2 + /?,,«, are assumed to be known, and v 1 ,v 2 are 
injection terms that are to be determined so that y — » y, v^v as time increases. 

The state estimation error dynamics is obtained by subtracting eq. (37) from eq. (36): 

\ee= e q+ { c \4?h - Cltfl ) “ V 1 


~ a 2\ e 6 + {~ a 2\ c l4 +a 24) 7 7l + a 23 e v +a 25 r l2 +c 25 7 /2 _v 2 


(38) 


\fv ~ a 3l e 0 + a 33 e v + ( _<2 31 c 14 +<2 34) 7 7l 

The injection terms are designed in a second order sliding mode control format, 7 in particular, in a super- 
twisting form 8 

. 1/2 


v 1 =a 1 \e 0 \ sign(e 0 ) + f\ J sign(e 0 ) 

i il/2 


v 2 = ^2 \ e a 


sign(e q ) + P 2 \sigri(e q ) 


(39) 


1 HI J 

e e = ye - ye > e q = y q - y q > *v= v -v 

The injection terms (39) guarantee convergence e 0 ,e q —>0 in finite time. Therefore, in the sliding mode we 
obtain 

c 147l— ^25*72 =Vl 

< (—^21^14 + a 24 ) r !l + a 23 e v + ^25% + ^25^2 = v 2 (40) 

/v = a 33 e v + ( _<2 31 c 14 + ^34 ) 7 7l + ^35% 

The bending mode estimates are derived from the error dynamics eq. (40) bearing in mind that . 

V 1 


The coefficient 


q 2 =- 

c 14 ~ c 25 

fl\ = [^23 {~ a 31 c l4 + a 34 ) + a 33 a 21 c 14 “ ^24^33 ] *7l + 

a 24 

~ — [^2 ~ a 33 v 2 +^25*72 ~( a 25 + a 33^2s) l l2 + a 33 a 25 r h\ 
A 24 

^77 = [ a 23 (“^31 c 14 + ^34 ) + a 33 a 21 c 14 “ ^24^33 ] 

a 24 


(41) 


(42) 


(43) 


in eq. (42) is positive, which means that eq. (42) is unstable. Instability of eq. (42) adds difficulty in finding a 
solution of this equation, which gives the estimate of the bending mode fj \ . In order to find a bounded solution of 
this unstable equation, the method of Stable System Center (SSC) 9 " 11 is employed. The main result of the method of 
SSC is formulated in the following theorem. 10 
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Theorem. Given unstable differential equation 


i=Qt+m, zj^R m (44) 

with the measurable input f(t ) , and the following set of conditions being met: 

(1) the matrix Q in eq. (19) is nonsingular and non-Hurwitz; 

(2) the input function f(t) can be piecewise represented by known linear exosystems with a characteristic 
polynomial given by 

P k (A) = A k + p k _\2 k 1 + ... + p x ) i + p 0 , (45) 

the state reference profile ^ e S J1 W can be generated by the matrix differential equation 

+c k J (k ~ 1) +... + cJ + c 0 4=-(p k _ l f +...+ PJ + PJ ) (46) 

where the constants c k _ 1 ,...,c l ,c 0 are chosen to provide any desired eigenvalues ( convergence rate) for ^ , 
and the matrices P k _ 1 ,...,P 1 ,P 0 e y\ mxm are given by 

P kA =[l+c kA Q 1 +... + c 0 2 k )[l + p k _ { Q 1 +... + PqQ k ) -/ 

P k - 2 = c k-iQi 1 +-"+ c o2i ik 0 ~{ p k-i+ b'iPk-iQi 1 +•••+ PoQi (k ") 

= c,Q~ l +c 0 Q~ 2 -{P k _ l+ iy( Pl Q- l +p 0 Q- 2 ) (47) 

P 0 =c 0 Q- l -{P k _ l+ l)-p 0 Q- 1 


The output filter (46) for the observer (42) is designed based on eqs. (45)-(47). 

x , 0.755.S 2 + 0.242s + 0.602 ,, 

ill (s) = 1 -45 f(s), 

s + 3 s + 3s + 1 

fit) = -3330v 2 - 44.666v 2 - 0.5/> 2 + 1 ,333/> 2 + 2.666 • 1 tr 4 // 2 

Next, the measurements (32) can be cleaned up using the estimated bending modes 
yOclean ~ y& ~ C \\P\ 


yqclean ~ y& ~ c 25 7 ll 

and the cleaned measurements are used in the control law (34) that becomes 


1 / •• 2 3 

~ \0 C +1 .75a) n 0 c ~ ^2\y Oclean — 1 -2 5co n y q C lean +2.15 0) n e^ + CO n 


where e 0 = 0 C - y 0clean . 


(48) 


(49) 


(50) 


E. Simulations of the compensated Ares I flight control system with the bending mode observer and cleaned 
up measurements 

The results of the simulations of the coupled flight control system (28)-(32) with the bending mode observer 
(37), (39), and (48), cleaned up measurements (49), the PID controller (50), and co n = 1 are shown in Figs. 7 and 8. 
The accuracy of the bending mode estimation is very high and is illustrated in Figs. 9, 10. 
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Fig. 7bZoomed output tracking in system with 
coupling, bending mode observer and cleaned 
measurements 
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Fig 8a The control function for system with coupling, exact 
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Fig 8b Zoomed control function for system with 
coupling, exact rate measurement and corrupted 
attitude measurement 


The output tracking with cleaned measurements (Fig. 9) is comparable to the one with exact measurements 
presented in Fig. 3. 
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Fig 9b Zoomed bending rate estimation 


Fig 9a The bending rate estimation 
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Fig 10a The bending position estimation 



Fig 10b Zoomed bending position estimation 


F. Discussion of the obtained results 
Advantages: 

1. Use of sliding mode observers improves significantly the robustness of flight control of Ares I. 

2. The proposed approach is capable of working in concert with any controller, improving its robustness. 

3. The bending mode dynamics are treated as the external disturbance while being estimated by a sliding mode 
disturbance observer. It makes the proposed algorithm robust to the parameters of the bending mode 
dynamics. 

4. The proposed technique allows controlling the stability margins. Assuming the sliding mode observer 
estimates the flexible modes (treated as the external disturbances) exactly, control system completely 
compensates for the flexible mode effects on sensed data and flexible coupling between the attitude control 
channels. Therefore, only the PID controller runs a feedback and defines the stability/phase margins. 

5. The proposed techniques can be easily implemented in a gain scheduling fashion. 

Disadvantages 

1. Using the proposed sliding mode disturbance observer, it is possible to estimate one bending mode and the 
bending mode rate only, since the number of estimated disturbances must be equal to the number of 
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measurements. The higher order bending modes and sloshing modes cannot be estimated using the sliding 
mode disturbance observer. 

2. The coupling that involves the sloshing modes cannot be compensated by means of the control function, since 
the sloshing modes cannot be observed using the proposed sliding mode disturbance observer. 

3. We cannot enjoy full robustness to the parameters of the bending mode dynamics while estimating bending 
modes unless the sliding mode disturbance observer compensated dynamics are minimum phase. In our case, 
the nonminimum phase dynamics required employment of the method of stable system center. This method 
requires partial knowledge of the bending mode dynamics. 

4. A possible disadvantage could be in high frequency control computational cycles to allow sliding mode 
observer to converge and provide a high accuracy of estimation. 

In the next section, another approach based on a full state sliding mode observer is proposed. This observer 
estimates the states of the Ares I dynamics altogether with the dynamics of bending and sloshing modes. 

The expected advantages of the newly proposed technique are in the possibilities 

(1) to address estimate multiple bending and sloshing modes; 

(2) not only to clean up the measurements but to fully compensate the coupling terms in CLV dynamics using 
feedforward patches in a simple PID controller. 

The possible disadvantage is in the requirement to know accurately the mathematical models of the CLV, bending, 
and sloshing modes. The deviation of the actual frequencies of bending and sloshing modes from the model values 
by up to 10% could be sufficient for obtaining reasonably accurate estimation and compensation. That yields the 
observer of a high order of that could be difficult to tune and is potentially sensitive to parametric uncertainty. 


VI. Flexible CLV attitude control using sliding mode state observer 

In this section we will design the state sliding mode observer and observer -based control while relaxing the 
assumptions A1 and A2. The new relaxed assumptions are 


A3. Multiple bending and sloshing modes are taken into account; 

A4. The dynamics of actuator are taken into account. 

The Ares I linearized pitch plane flight dynamics are defined by eqs. (25): 
jx = Ax + BU 

1 y = Cx 

and the actuator dynamics are presented in a form of LTI system of differential equation (26): 

Jz = Qjz + Q 2 x + (3u 

jt/ = Q 3 z 


(51) 


(52) 


where 


x = 


<9, q , v, //, , m , % - 7 h > 8 \ > s 2 > $3 - S A 


r T 

x 2 


□ 11 is a state vector of the Ares I pitch plane dynamics, 


z e □ 5 is a state vector of the actuator, 

U e □ 6 is the output of the actuator that is fed to the Ares I dynamics in eq. (51) through the distribution matrix 

B g □ 11x6 , and u efl 1 is the control function that is the servo command input of the actuator (52). 

Remark 4 . Since the entries of the vector U represent the outputs of the actuator (some of them may be equal to 
zeros) that are acting on the inputs of the CLV, the distribution matrix B in eq. (51) consists of “1” or “0” elements 
only. 

The distribution matrix B is defined as follows: 
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B = 


0 0 
1 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
1 0 0 
0 0 0 
0 1 0 
0 0 0 
0 0 1 
0 0 0 
0 0 0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


A. PID controller design 

The PID controller is designed for the simplified flexible CLV attitude dynamics given by eq. (25) that takes 
into account two bending and two sloshing modes. The tracking error dynamics are obtained after 2 differentiations 
of the ideal tracking error e q =6 c -6 while neglecting the actuator dynamics and assuming perfect measurements 
(this assumption is not realistic, since the measurements will be cleaned up). Therefore, 

u = ——(o c +lJ5co n d c -U5co n q + 2.l5o}^eg+ 0 }llegdt- x ¥(d,v,Tj l ,...,rj 4 ,Si,...,S 4 )\ (53) 

b 2 1 V ’ 

where Z?21 * s taken as in eq. (29), and 

¥(<9,^, 7ft, 772,773, 774, S h S 2 , S 3 , S 4 ) = 

a 2,l^ + a 2,2 ( l + a 2,2 v + a 2A r ll +a 2,5 r l2 +a 2,6% + a 21 r l4 + a 2,^\ + a 2,9^2 + ^ 2 , 10^3 + ^ 2 , 1 1^4 
with the coefficients a 2 j j = 1, 2, ..., 1 1 taken from the second row of the matrix A . 

Remark 5. It is worth noting that compensating of the coupling terms collected in the function (54) becomes 
possible due to observing of all states of system (25) via a sliding mode observer that is presented later. 


B. Coordinate transformation 

The mathematical model of a flexible CLV in eq. (51) is rewritten in a form suitable for the sliding mode 
observer design, where the first part of the state vector coincides with available measurements, i.e. with vector y : 

j*l = A n x l +A n x 1 +B l U 

[^2 — Ar£\X\ + A22X2 + B2JJ 

y = C 1 X 1 + C 2 x 2 —> Xi = Q 1 (y — C 2 x 2 ) , C\ = I 2x 2 (55b) 

where A n , A 12 , A 21 , A 22 , £] , # 2 , C] , C 2 are P ar tial blocks of original matrices of eq. (51) of corresponding 
dimension. 


The goal of the basis transformation is to present system's dynamics in new basis (y,x 2 ) T as: 

\y - \\y+ \2 X 2 + ^i u 
\*2 = ^21 y + ^22 x 2 + B2U 

where 


(56) 
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(57) 


A x , = (C 2 A 21 + C ] A ] j ) Q 1 

A 2 = C 2^22 + C lA2 _ ( C 2^21 + C lAl ) C 1 ' C 2 

A 2 i = A 21 Q 1 

a 22 = a 22 - a 21 q 1 c 2 

+ C 2 B 2 

Z?2 — /^2 

Remark 6. In order to reduce complexity of the calculations while designing the sliding mode observer for system 
(56) we take into account the actuator dynamics (52) in a special way. The designed scalar control function u in 
eqs. (53), (54) is assumed known and is to be fed through the actuator to generate the vector control U . This control 
U is, in fact, the vector-output of the actuator and also assumed known. Assuming that actuator output is known 
means assuming there are sensors on the vehicle to measure them. This can add complexity to the controller 
implementation. The vector control U is then distributed among flexible CLV dynamic equations (55a) by means 
of known matrices B^B 2 , and the products B^U , B 2 U are assumed known while designing the state sliding mode 
observer. 


C. State sliding mode observer design 

The state sliding mode observer is designed based on the Ares-I dynamics written in a form (31). Introduce the 
observer's dynamics: 

[y = A n y + A n k 2 +B x U + g 

[*2 = A 2l9 + A 22 x 2 +B 2 U + Lg 
where y,x 2 are the observer’s states. 

The design matrix L e □ 9x2 is selected to ensure the desired rate of convergence e y , e x ^ — » 0 in the sliding 

mode. The sliding mode injection term g = [g 1? g 2 ] 2 G 0 2 that drives the estimation errors 
e y = y~y> e x 2 = x 2 ~*2 to zero (e y , e x ^ — » 0 ) as time increases is designed in accordance with super- twisting 
algorithm. 7,8 The vector-control input U and the output y are available for measurements. 

The observation error dynamics are identified 


e y - \ l e y + \ 2 e x z - 8 ( e y) 


x 2 


^2l e y + A 22 e x 2 -L-g(e y ) 


(59) 


The injection super-twisting terms are designed 

il / 2 


g, (e y] ) = a x |e yl | sign (e y] ) + /?,{ sign(e y] ) d r 
S2i e yl) = «2 k^2| 1/2 *ign[e y2 } + P 2 \sign(e y2 }dz 


(60) 


Remark 7. The super-twisting injection terms (60) are continuous functions. 


The observation error dynamics in the sliding mode become 


J 0 = A2% -g< e y) 

\% =^22 e x 2 ~L-g(e y ) 


( A 22 e x 2 ~ 7A| 2 ) e x ^ 


(61) 


The design matrix L is to be selected to provide placement of desired eigenvalues of the matrix A 22 e x ^ - L\ 2 in 

eq. (51). For the purpose of simulation, matrix L was evaluated to provide a pole distribution according to the 
Butter worth polynomial with co = 5 . 

Finally, as soon as x 2 converges to x 2 , the full state vector X can be reconstructed as 
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iJcrVcA)]. (62) 

x 2 

Thereafter, the estimated values of the states x in (52) are to be used in the control law (53), (54). 

Remark 8. The estimated state sub-vector jq = Cf 1 ( y-C 2 x 2 ) = y-C 2 x 2 is, in fact, the cleaned up measurement 
vector y . 


D. PID controller augmented by state sliding mode observer. 

In summary, the linearized pitch plane dynamics of a flexible CLV that takes into account 2 bending and 2 
sloshing modes as well as the actuator dynamics given by eqs. (51), (52) 
f x = Ax + BU fz = Q.iZ + Q 2 x + f3u 

{ y = Cx | U = Cl 3 z 

are controlled by the control law 

u = -^—(& c +1.15a) n d c - U5co n q-\-2.l5co^e 0 + cD^\e^dt- x ¥(6,v,q l ,...,q 4 ,S l ,...,S 4 )\ (63) 

*21 V ’ 


'V(0,q, v, a ;, , f] 2 , f] 3 , tj 4 , ,S 2 ,S 3 ,S 4 ) = 

a 2j@ + a 2,20 + a 2,3^ + <2 2,4 / 7l + a 2,502 +a 2.6 ? h + a 27*?4 + a 2,8^l + a 2,9$2 + a 2,IO*3 +a 2A 1*4 


with the states x= <9, q , v, rji , rj 2 , rj 2 , 774 , ^ , S 2 , S 2 , S 4 estimated by the state sliding mode observer. 


E. Simulations of the compensated flight control system of Ares I with two bending modes, two sloshing 
modes, actuator dynamics, and the sliding mode observer. 

The results of the of the simulations of the compensated flight control system (51), (52), (53)-(64) that 
comprises the Ares I pitch plane dynamics linearized at time t = 60 sec , two bending modes, two sloshing modes, 
actuator dynamics, and the sliding mode observer are presented in Figures 11-12. The following non-zero initial 
conditions are assigned to the state variables during the simulations. Angle states are given in radians, velocity is in 
feet per second, bending states are in units of generalized (modal) displacement, and slosh displacements are in feet. 

x = (0,q,v,q l ,r/ 2 ,r] 3 ,ri 4 ,S l ,S2,S 3 ,S 4 ) T = 

T (65) 

(lo^ao^io^-io^s-itrVio^s-KT^io^-icrVio^-io" 2 ) 

The time history of estimating the first bending mode and first bending mode rate are shown in Figures 1 1 and 
12. The time history of estimating the second bending mode and second bending mode rate are shown in Figures 13 
and 14. The time history of estimating the first sloshing mode and first sloshing mode rate are shown in Figures 15 
and 16. The time history of estimating the second sloshing mode and second sloshing mode rate are shown in 
Figures 17 and 18. After short transient responses ~ 4-5sec the estimation accuracy is almost perfect. The time 
history of measured and “cleaned up” (estimated) pitch angle 0 and pitch rate q are shown in Figures 19 and 20. 
After a short transient response ~ l-2sec the measured values, corrupted by the bending modes, are accurately 
cleaned up. The time history of estimating the translational velocity is shown in Figure 21. After a short transient 
response ~ 4 - 5 sec the estimation accuracy is almost perfect. The pitch angle tracking and a corresponding control 
function are shown in Figures 22 and 23. Practically no effects of bending and sloshing modes are observed. 

In order to study effects of the external disturbances the flight control system was simulated upon the 

_o 

instantaneous change of the initial conditions for the pitch angle by -3-10 at t = 40 sec . This initial condition 
approximately represents the vehicle state immediately after an external disturbance or attitude command such as a 
wind gust or maneuvering. The results of the simulations are shown in Figure 24. The pitch angle tracking recovers 
after a short transient ( ~ 4 - 5 sec ). 
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VII. Conclusions 

The attitude control dynamics of a flexible launch vehicle model are studied while taking into account bending 
and sloshing modes. The linearized mathematical model of the flexible Ares I CLV was presented, which took into 
account two bending and two slosh modes as well as the actuator dynamics. 

A simple PID controller was designed and simulated for the vehicle dynamics model linearized at 60 sec of 
flight assuming (a) no coupling between the vehicle attitude dynamics and translational motion; (b) both two slosh 
and bending modes as well as actuator dynamics were neglected; (c) effects of the bending modes on the attitude 
sensors were neglected. 

First, the compensated vehicle dynamics with a simple PID controller was studied for the Ares I model that 
included two bending modes, two slosh modes and the actuator dynamics, without compensation for the bending 
dynamics. The resultant behavior of the Ares I compensated by the standalone PID controller appeared to be 
unacceptable due to classic control- structure feedback instability. 

In order to improve the flexible CLV attitude control system robustness to flexible body bending dynamic and 
propellant slosh modes, the sliding mode observers were designed. 

The first sliding mode observer treated bending and sloshing modes as external disturbances. Only first bending 
mode and its rate were possible to estimate using this type of observer. As soon as the bending modes are estimated 
their effects on the attitude and attitude rate measurements were compensated. The cleaned up measurements were 
used by the controller. The coupling terms that describe influence of bending and sloshing modes on the Ares I 
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dynamics were not compensated due to a lack of necessary estimates. However, it was demonstrated that the Ares I 
dynamics could be stabilized without using notch or phasing filters. 

The second sliding mode observer proposed in this research is capable of estimating bending and sloshing 
modes together with the states of the flexible CLV. Effects of two flexible modes and two slosh modes on 
measurements and the coupling effects have been studied and compensated. The intensive simulation study 
demonstrated stable, accurate performance of the flight control system with the augmented PID controller. The 
compensated vehicle dynamics again became stable without using notch or phasing filters. The Ares I compensated 
dynamics demonstrates practically no bending and sloshing influence and is comparable to the rigid body 
compensated dynamics. 

The proposed observers can operate with any controller and do not rely on classical bending filters. The 
application of the proposed observation technique can be accomplished in a gain- scheduled fashion to accommodate 
time- varying plant parameters. 

Our approach assumed a linear perturbation model of the CLV, which effectively represents the actual behavior 
of the present vehicle under most operating conditions. 19 In order to address nonlinear phenomena of the flexible 
CLV attitude dynamics, a nonlinear model could be used for designing the observers. Corresponding higher order 
sliding mode observation techniques have been developed recently. 12,18 
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